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Abstract. We present a search for neutral supersymmetric Higgs bosons decaying to t + t~ pairs 
produced in pp collisions at y/s = 1.96 TeV. The data have been collected with the CDF II and DO 
detectors at the Tevatron collider at Fermilab (1 fb _1 of integrated luminosity per experiment). No 
significant excess above the standard model backgrounds is observed. We set exclusion limits on 
the Higgs production cross-section times the branching fraction of its decay to t + t~ pairs for Higgs 
masses in the range from 90 to 250 GeV/c 2 . We also set exclusion limits on MSSM parameters vtla 
and tan (3 in several benchmark scenarios. 

PACS. 14.80.Cp Non-standard-model Higgs bosons - 13.85.Rm Limits on production of particles 
- 12.60.Fr Extensions of electroweak Higgs sector - 12.60.Jv Supersymmetric models 



1 Introduction 

The Higgs mechanism provides the scheme for the elec- 
troweak symmetry breaking, that is accepted in the 
standard model (SM) as well as some of its supersym- 
metric extensions, such as Minimal Supersymmetric 
Standard Model (MSSM) [lj. The Higgs boson, if it 
exists, remains to be found, and the Tevatron exper- 
iments are actively searching for it. Due to low pro- 
duction cross-section predicted for Higgs in the SM, 
the SM Higgs cannot be observed with the amount of 
data collected to date. The production of MSSM Higgs 
may be enhanced relative to that of the SM by several 
orders of magnitude for certain areas of MSSM param- 
eter space [2], making discovery of a MSSM Higgs at 
the Tevatron possible with less integrated luminosity. 

In MSSM, the Higgs sector is comprised of two dou- 
blets of complex scalar fields. These correspond to the 
five physical particles: two charged bosons i^ ± , one 
neutral CP-odd boson A and two more neutral CP- 
even h (light) and H (heavy) particles. At the tree 
level, two free parameters are sufficient to describe 
properties of these bosons. Most often, these param- 
eters are chosen to be the mass of the CP-odd state 
rriA, and the ratio of vacuum expectation values of 
Higgs coupling to down-type and to up-type fermions 
denoted tan (3. At low tua and high tan /3 the A is al- 
most mass-degenerate with either h or H , and the pro- 
duction cross-section of the <j) is enhanced by tan 2 /3. 
It is conventional to use <ft to denote any of the A, h 
or H. 

The MSSM Higgs at the Tevatron is expected to 
be primarily produced in the gluon fusion process, and 
the second most significant process is bb — > (j). It de- 
cays roughly 90% (10%) of times into bb(rr) pair. The 
analyses presented in this paper are devoted to the lat- 
ter channel. CDF has searched for the Higgs decaying 
into the final states r e Th a d, T^Thad and r e r M , where the 



subscript denotes either the corresponding leptonic, or 
a hadronic r decay. The DO analysis uses the r^Thad 
final state. 



2 Data and Monte Carlo Samples 

The analyses are based on the data collected at the 
Tevatron in pp collisions at \/s = 1.96 TeV. Both 
CDF and DO experiments used data samples that cor- 
respond to the integrated luminosity of ~1 fb _1 . 

The data of CDF have been collected using di- 
lepton triggers for the T e T^ final state, and lepton+ track 
triggers. At DO, single muon triggers have been used, 
that require hits in the muon system in conjunction 
with a high-momentum track in the central tracking 
volume. 

For background estimates and efficiency calcula- 
tion both experiments have used Monte Carlo samples 
created with PYTHIA event generator [3] and GEANT 
detector simulation [3]. 



3 Backgrounds 

The largest and irreducible background in this search 
comes from events with Z/^y* — > rr. This background 
is estimated using Monte Carlo simulation. 

The QCD multi-jet events and W+jet events con- 
stitute the second largest source of the background. 
These contributions can be controlled with well cho- 
sen selection criteria. 

Finally, there is a number of small backgrounds 
that are taken into account, such as di-boson and tt 
production, and Z — > ee, jifi. 
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Table 1. Event count in data and predicted background 
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Data 
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1000 


374 


1144 



4 Analysis 

4.1 Event Reconstruction and Selection 

The analysis of each experiment is only briefly de- 
scribed in this paper, a complete report can be found 
in [5] and [6] for the CDF and DO, respectively. 

For the reconstruction of a leptonic tau decay trig- 
ger leptons are used both at CDF and DO. These lep- 
tons are required to match a high-momentum track. 
They also have to pass isolation requirements on the 
total amount of energy deposited in the calorimeter 
and scalar sum of charged track momenta detected by 
the tracking system in a cone around the leptons. 

The r leptons decaying hadronically are identified 
by detection of collimated jets. The primary r decay 
channels contributing to this signature are the 7r ± ^ r , 
7r ± 7r°^ r and 7r ± 7r^7r ± (n7r°)^ r . The jets detected in 
calorimeter are thus matched to one or three charged 
tracks, the angular size of the jets is limited to a nar- 
row cone and, in case of CDF, isolation requirements 
on detected particles outside of the r jet cone are im- 
posed. The DO experiment identifies r leptons with 
neural networks developed on Z — > rr events. With 
such selection, multi-jet QCD backgrounds are sup- 
pressed. 

Other backgrounds are reduced by several require- 
ments at the event level. Events at DO have to pass 
W-veto that removes W — > [w events, with mw re- 
constructed using the momentum of the muon and 
the vector of the missing energy. Similar backgrounds 
are suppressed by CDF with the requirement that re- 
jects events where the missing energy points away from 
the combined direction of the reconstructed rr decay 
products. Other methods used by CDF and DO to elim- 
inate lesser backgrounds can be found in [5] and [6]. 

After applying event reconstruction and selection 
to the data samples, both CDF and DO find that the 
observed yields in the data match the expected counts 
of background events within uncertainties (see Table 
[1]). Note that CDF errors are statistical only, while the 
errors of DO also include systematics. 



4.2 Results for a(pp -> 0X) Br(0 -> rr) 

The selected data samples have been searched for the 
signs of the Higgs signal in the mass range from 90 to 
250(200) GeV/c 2 by CDF(DO). 

CDF used the binned likelihood fitting technique 
looking for a peak in the distribution of the invariant 
mass of the visible rr decay products and the vector 
of the missing energy. A small excess of events has 
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Fig. 1. The limit on cr{pp <j>X) • Br{4> rr) from CDF 
as a function of Higgs boson mass. 
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Fig. 2. The limit on a(pp -> (j)X) • Br((j) -> rr) from DO 
as a function of Higgs boson mass. 



been observed in the r e r M sample at the mass of about 
140 GeV/c 2 . The statistical significance of the excess 
is approximately two standard deviations. 

Neural network analysis has been developed by DO 
and used event parameters such as the visible mass 
(defined as at CDF), transverse momenta and pseu- 
dorapidities of the decay products. The distribution of 
the neutral network output for the data has been com- 
pared to the expectation for the background, and no 
statistically significant deviation has been observed. 

The data of both CDF and DO are consistent with 
background-only observation. Consequently, both ex- 
periments set a 95% C.L. exclusion limit on a(pp — > 
cj)X)'Br((j) — > rr), the product of the Higgs boson pro- 
duction cross-section and the decay rate of the Higgs 
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|j = -200GeV, M 2 = 200 GeV, m. = 0.8 M SUSY 
M SUSY = 1 TeV, X t = V6 M SUSY (m h max ), X t = (no-mixing) 
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Fig. 3. CDF exclusion area for tan f3 as a function of tua for 



M = +200 GeV, M 2 = 200 GeV, m g = 0.8 M SUSY 
M SUSY = 1 TeV, X t = V6 M SUSY (m h max ); M SUSY = 2 TeV, X t = (no-mixing) 
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several benchmark scenarios overlaid with the LEP results 



to r-pair. This limit is dependent on the Higgs boson 
mass. The limit from CDF is shown in Fig. dJ The ob- 
served exclusion is slightly weaker than expected for 
given statistics due to the 2a event excess mentioned 
above. The limit from DO is presented in Fig. O 



5 Interpretation of Results within MSSM 



The measured limits on the a(pp — > <j>X) • Br((j) — > rr) 
can be used to derive exclusions on the free parame- 
ters of the MSSM. As the number of the parameters 
is large, the common approach to interpreting results 
of Higgs searches is to fix most of the MSSM param- 
eters to agreed-upon values and then derive the ex- 
clusion regions in the plane (ra^, tan/3), the two pa- 
rameters most directly related to MSSM Higgs prop- 
erties. There are several ways to fix MSSM parameters 
corresponding to the most indicative cases. These are 
called MSSM benchmark scenarios [7J. Out of the sce- 
narios described in [7] , CDF and DO employ the m™ ax 
and the no-mixing scenario with the case of the mix- 
ing parameter of the Higgs doublets \i being positive 
and negative. The m™ ax scenario is defined as the set 
of MSSM parameters that maximizes the value of the 
mass of the lightest Higgs boson. The second scenario, 
no-mixing, is similar to m™ ax , but with the mixing 
parameter for the stop quarks set to zero. Total of four 
cases are considered. 

The CDF exclusion regions for the four scenar- 
ios are drawn in Fig. [3] and the results from DO are 
found in Fig. [H and O In the studied mass region the 
tan/3 values above 40-60 are excluded in the consid- 
ered MSSM scenarios. While no signal is observed and 
the data are consistent with the SM backgrounds, the 
measurements presented here provide the most con- 
straining to date limits on Higgs decays to rr channel. 
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Fig. 4. DO exclusion area for tan /3 as a function of for benchmark scenarios with fi < overlaid with the LEP 
results and previous MSSM Higgs searches at the Tevatron with lower luminosity 
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Fig. 5. DO exclusion area for tan (3 as a function of for benchmark scenarios with \i > overlaid with the LEP 
results and previous MSSM Higgs searches at the Tevatron with lower luminosity 



